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El problema del análisis de la flexión simple en piezas rectas de sección constante con ciertas 
hipótesis complementarias y otras condiciones de forma y solicitación nos conduce a un modelo regido por 
una ecuación diferencial lineal de coeficientes constantes con condiciones de contorno que puede ser 
reducido a un problema unidimensional definido en un intervalo. Este problema es, como se sabe, 
fundamental en el campo de la edificación por el predominio de estas piezas en arquitectura e ingeniería. 
La presencia de cargas y momentos puntuales en las solicitaciones de vigas es imprescindible para 
completar un modelo operativo ya que permite ampliar substancialmente el campo del análisis a vigas de 
varios tramos y otros problemas de interés en el contexto estructural. Pero esta presencia nos impide 
efectuar holgadamente integraciones "convencionales" con lo que se necesita acudir a otros procedimientos 
matemáticos entre los que se encuentran los métodos numéricos. 

El método de integración por diferencias finitas es un método numérico universal y 
extraordinariamente eficaz que resulta una solución idónea por su simplicidad y adaptabilidad a cualquier 
tipo de planteamientos en este ámbito. El sistema de ecuaciones al que conduce la resolución de los 
problemas puede hoy día ser resuelto con gran facilidad por la mayor parte de calculadoras ordinarias y 
programas informáticos. 

Después de una exposición del procedimiento con carácter general descrito sobre los casos más 
simples para ecuaciones de 2” y 4” orden, se ha procedido a configurar unas fórmulas de uso sistemático 
totalmente adaptadas a las condiciones del problema de la flexión simple en lo que hemos denominado 
"métodos en dos etapas” partiendo de la caracterización de las condiciones de conjugación de la unión 
rígida, el tratamiento de las secciones extremas y, en su caso, de la unión articulada. Con esta descripción 
pormenorizada puede ya procederse reiteradamente sobre las piezas de uno o varios tramos sean cuales 
sean sus vinculaciones de extremo. Se ha dado la mayor importancia a los ejercicios prácticos que son los 
que permiten comprender la metodología operativa y obtener los resultados buscados en los distintos casos 
comparando sus soluciones con las obtenidas por otros métodos. 

Finalmente somos conscientes de que en un trabajo de estas características es inevitable la aparición de 
erratas o de falta de precisión que iremos corrigiendo en revisiones posteriores. 


Pedro Galán del Sastre 
Ramón J. Zoido 
Mayo 2014 
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1 Las Diferencias Finitas. 
Aproximación de las derivadas por 
diferencias 


1.1. Diferencias finitas entre los valores de una 
función. Operadores 


Una forma de discretizar de manera uniforme una variable real continua x es considerar 
los valores 1, = tp + 1h, con 1 € Z, a partir de un cierto valor xy. De esta manera definimos 
una red unidimensional uniforme de paso h constante, R|, formada por los puntos 2, de la 
variable x, es decir, 

Rs, = [Lil iez 

Si consideramos ahora la función real de variable real, continua, f(x) en un intervalo 
J = [0,1] CR y la variable x discretizada de esta manera, la función restringida en J a 
los valores sobre esta retícula discreta de puntos R;, [0, l] determina una sucesión numérica 


Ufibocien" 


LSibocien —= A A 


con = Ple), du =%h para 0<3< Ny N= _: que es lo que podemos denominar 
proyección de f (x) sobre la red R! [0,1] (figura 1.1). 

Así, la función f (1) = 1% — x definida y continua en todo R, fijado 27 = 0 y paso h = 1, 
proyecta la sucesión ([f;),.7 =41..., —6,0, fo = 0,0,6,24, 60,... ). 


Figura 1.1: Representación discreta de una función. 


1 Las Diferencias Finitas. Aproximación de las derivadas por diferencias 


Definimos ahora, para todo 1 € Z, 


Anti = Fira — fa, (1.1) 
que se lee diferencia primera progresiva de f;. 
Puede formarse otra sucesión: 


ARA EA A A A 


y, en ella puede ser reiterada la definición (1.1) obteniendo: 


Ap(Apfi) = Apfi = Aplfiva — 4) = Apfira — Apt = Fira — 2fi41 + fi 
o diferencia segunda de f; en sentido progresivo. Así puede formarse la sucesión de diferencias 
segundas: 


ds LEL , Añfza, Añto, Añf1, Año, a A ... ) 


y, en general la diferencia de orden k: 
k k po 
Fi+r — ( 1 ) Fekca + ( 2 ) uh A [=P ( j ) Fi+kj Aa (1 fi 


que determinará la sucesión da pe DES AT AB, AE fa, 2. AB, Aria nd » 

En determinadas aplicaciones aparecen problemas regidos por ecuaciones en diferencias 
finitas de orden n con condiciones iniciales y/o de contorno. Las ecuaciones lineales con coefi- 
cientes constantes tienen un tratamiento similar al de las ecuaciones diferenciales ordinarias. 

En (1.1) queda implícitamente definido el operador A, como un operador lineal llamado 
operador en diferencias progresivas. Análogamente podemos definir las diferencias regresivas 
en la sucesión «f;) ¿7 por 


Arfi= fin fia, 


y, asimismo establecer las diferencias medias o centrales: 


Amfi= 3 (Apfi+ Arf) => fr fica). 


Pueden definirse otros operadores sobre la sucesión + f;) ¡ez: Piendo 1 el operador identidad, 
If; = fi; el de traslación T se define por T'f, = f,,1 y su inverso T"*f, = f;_,; ambos nos 
permiten ya establecer algunas relaciones entre ellos: 


Artio = fi fi=T ff =(1-0D fi, 
Artio = fifa=f-T*f= (1 =T”>) Fa, 
por lo que podemos escribir, por ejemplo, Ay = T — T o bien A, =1—T”?, 
Definidas la suma y la composición de operadores 
(01+02)f = Of+052f, 
(O, 0.02) f 01 (02f), 


1.1 Diferencias finitas entre los valores de una función. Operadores 


y definido el producto escalar de un número A € R por un operador: 


MO)f=A(0f), 


resultan inmediatas las siguientes relaciones entre los operadores T, T'7*!, Dig 


Ap=T -—I, T=eIFAs 
Aj=I=T", PA=T= An 
Tobr=T =1=D PAGA SÓ 


Estos operadores satisfacen las leyes conmutativa y asociativa ordinarias para la suma 
y el producto siendo, además, el producto distributivo respecto de la suma. Sin entrar en 
su estructura y utilidad baste decir que careciendo de sentido independiente de la función 
que debe de ser colocada a su derecha, pueden, no obstante, ser tratados como cantidades 
algebraicas y formar polinomios que constituyen, a su vez, operadores más complejos. En 
nuestro contexto, estos operadores pueden ser relacionados con el operador derivada: 


D"f(c)=£" (2) 


para funcionas analíticas (de variable real) f (x) a través de la fórmula de Taylor: 


Tf()=f(a+h) 
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por lo que: 
I+A,=T=e'P, 


y por lo tanto: 


D= Ino T= In (14 Ay). 

Podemos extender los procedimientos de proyección sobre redes uniformes discretas a las 
funciones f (1) segmentariamente continuas en el intervalo J = [0,/] € R, con un número 
finito de discontinuidades de primera especie o de “salto finito” en ./. Por el momento, en los 
nodos x = zx; de la proyección de f (x) que coincidan con discontinuidades de la función con- 
vendremos en que eventualmente la función adopta el valor promedio o, en general, podemos 


considerar por definición que la proyección de la función sobre la red es de la forma: 


Pu) =f=3 0 (0) +f(07)), 


N]|r 


siendo 1, y x+ el límite de la función f en x; por la izquierda y derecha respectivamente, 


lo que equivale a substituirla por una función continua mejor o peor aproximada en función 
del paso h de R;,. 
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1.2. Aproximación de derivadas por diferencias finitas 


Para una discretización ER los valores de la variable x: x, € R¿[0, 1], con 2, = 0 + 1h, 
¡eK=ZN[0,N] y N=!, se trata ahora de aproximar yW (x;), k=1,2,3,...,n, de una 
función y (1) por dterendías entre los valores de [Yid:ex = LY (Ei) herx conociendo cualita- 
tivamente el error que dicha aproximación comporta y que será el error de discretización. 
Este error será normalmente puesto de manifiesto por una comparación con el desarrollo en 
serie de Taylor de la función en los entornos de x, por lo que, en lo que sigue, convenimos 
en suponer la suficiente derivación sucesiva de la función y(x) así como la continuidad de 
estas derivadas. Para poner de manifiesto este error se puede utilizar el llamado término 
complementario de Lagrange. Otras veces se señalará simplemente el orden de aproximación 
con la notación O (A”). 

En general la notación O (1”) denota una función en la variable x, así, cuando se escribe 
que la función f(x) = O (x”), por definición, la función f verifica que 


== GER, 


Es muy común utilizar el término complementario de Lagrange junto con el polinomio de 
Taylor para expresar cualquier función analítica como 


alo da PU A a o 
dd 3 hOA + a ÓS (n +1)! pa 


flo+h) = f(2)+f (2)h 
= Ta(f123h)+Ra ee h), 
donde 0 € (0, 1). Alternativamente, podemos expresar 
fl+h) =T,(f,23h) +0 (a+) 


puesto que 
16 En (L 0h) _ LD (0) 
E (n+1)! 

Las aproximaciones en diferencias de las derivadas no deben de ser vistas fuera del con- 
texto del problema numérico específico a resolver puesto que cada problema considerará un 
error admisible del que no sólo forma parte el de discretización. Con vistas a aproximar de 
una forma general soluciones de ecuaciones diferenciales lineales con condiciones de contorno, 
dispondremos aquí básicamente de aproximaciones en diferencias centrales de las correspon- 
dientes derivadas con órdenes de aproximación O (h?). Posteriormente realizaremos algunas 
correcciones propias del problema que nos A 


Si en la expresión y' = y' (1) = Xím a y para un valor Ax = h substituimos y, = y” (1;) 
ApYi A-Yi 
or o bien por obtendremos 
PO “Ag PO “Ag” 
/ Yi41 — Yi 
y: = y (2;) 7, 
o bien 
/ Yi — Yi-1 
== L¡) = y 
y; = y (25) » 


1.2 Aproximación de derivadas por diferencias finitas 


consecuencia de haber aproximado, para un cierto paso de red h fijo, la derivada por el 
cociente finito con diferencias progresivas o regresivas respectivamente. Ambas son aproxi- 
maciones O (h) ya que escribiendo el desarrollo de Taylor en el entorno de x;: 


h? h3 hi. R$ 
Yi = Ya hy Y gr pu + py +0 (A) (1.2) 
de aquí: 
+1 — Yi h 0 AS a hi 
mo =y; + 30 He Y + + EY +0 (1?) (1.3) 
así que 
y = ELE 40 (h). 
Así mismo, se puede escribir: 
h? AR O 
== a 01) (1.4) 
de donde: » a , 
Yi Yi-1 h " h mM h iv h v 5 
así que 


y _ Yi > +0 (h) 


Se mejora substancialmente la aproximación hasta O (h*) sumando (1.3) y (1.5), obtenien- 


do 
VYi41 — Yi me RE do 
or L= y + O (no) (1.6) 


así que 
rn Yi41— Yi-1 2 
'= KA F O (AS). 
Y; an + (17) 


Esta aproximación por cociente por diferencias medias o centrales, va a constituir nuestra 
aproximación tipo O (h?) de la derivada primera. 
Si en vez de sumar, restamos (1.3) y (1.5) obtenemos 


3 


h3. 
= hy! + 2 +0 (A?) 


Yi+1 — 2Yi + Yi-1 
h 


así que: 
a a 
y _ Yizi == Yi-1 PO (1?) 
obteniendo la aproximación O (h?) para y/. 
Para obtener una aproximación del mismo orden en la siguiente derivada deberemos con- 
siderar los desarrollos: 


ae ) 20 mM 25h" 1U 2h 0) Pa 
2 die 31 YT Al Y, 51 Y, 


Viso = Y +F2hy, A 5 Yy+O (7) (1.7) 
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y 
92)? 93)3 94p1 mn 95p? ” 9678 ce 
Yi-2 = Yi — 2hy, + 5) Ue == ar e A Yi — a Y ES MTS 0 (m7) (1.8) 
y restando ambas expresiones 
2eR> 2%hP 
Yi+2 — Yi-2 = 4hy! +2 al y +2 A O (h*) (1.9) 
De (1.6) podemos obtener: 
/ he YY h? VU Yé 
Yi+1 — Yi-1 = 2hy, 4 25 y; 7 25 y +0 (hn) (1.10) 


y, restando finalmente de (1.9) la ecuación (1.10) multiplicada por dos, obtenemos: 
'l 
Yi+a — Yi-2 — 2 (Yi+1 — Yi-1) = 2h y," + SUN +0 (17) 


que puede escribirse como 


Vara = ZYiga + Yi = Yi 
2h3 


=y + UN +0 (h*) 


así que 


y" = Yi+2 — 2Yig1 + 2Yi-1 — Yi-2 
! 2h3 
obteniendo la aproximación O (Ah?) para la tercera derivada. 
Sumando ahora (1.7) y (1.8): 


HO (AS), 


2 2 
Yi+2 + Yao = 2yi + 4h*y, +2 q A AA O (h3) (1.11) 
por otra parte la suma de (1.2) y (1.4) es: 
2,11 h? tu no vi 8 
Ya + Y 1 = 2, +hR%y, + 2% + 2 +0 (h ) (1.12 


para eliminar y/ entre ambas expresiones basta ahora restar de (1.11) la ecuación (1.12) 
multiplicada por cuatro. Se obtiene: 


: 1 E 
Yi+2 + Yi-a — 4 (Yir1 + Yi-1) = —6y, + Hoyo + Y +0 (m3) 


y de aquí: 
2 us 4 2 + 6 a 4 i— + 1— 1U il vi 
Yi+2 Yi+1 A Yi-1 T Yi-2 = y + y) +0(h%) 


por lo que podemos escribir 


su Via — Ús + Oye — ia E Y a 
Y; — he 


+0(h?) , 


obteniendo la aproximación O (h?) para la derivada cuarta. 


1.3 Resolución numérica de ecuaciones diferenciales ordinarias lineales 


1.3. Resolución numérica de ecuaciones diferenciales 
ordinarias lineales con condiciones en los extremos 
por aproximación en diferencias finitas 


Aunque no supone una restricción al método, nos referiremos a ecuaciones lineales de 4% 
orden con coeficientes constantes, es decir, de la forma: 


ay" (1) + azy” (2) + a2y” (1) + ary! (2) + a0y (1) = f (2), (1.13) 


con aj ER, f € C* ([0, 1]) y condiciones suficientes en los extremos del intervalo [0,1] para 
definir en él la solución particular buscada. 

Elegido un cierto número natural NV, se discretizan los valores de la variable x obteniendo 
la red unidimensional de paso h = » representada por el conjunto de puntos de RI [0,1]: 
zx, =0+1h,1€ ZNM|0, N], todos ellos pertenecientes al intervalo [0, 1]. 

En lo que sigue, a los valores que obtendremos mediante la correspondiente aproximación 
en los puntos x; los escribiremos mediante ys, Y, Y, y!” e y”, es decir, Y = Y (24), y, = y (24), 
etc... (a diferencia de la $ 1.2, donde esta notación la utilizábamos para los valores exactos 
de la función y, así como sus derivadas, en los nodos 2;) 


Se procede entonces a substituir en la ecuación de campo (1.13) las derivadas de los 
diferentes órdenes por sus aproximaciones en términos de diferencias entre los valores de 
R] [0,1]. Utilizaremos las aproximaciones centrales O (h?) obtenidas anteriormente: 


id e PA 
j 2h 
e io Zo Us 
Yi = 72 
mo — Vi+2 — 2Yi Ed 2Yi 1— Yi;i-2 (1.14) 
Y; 23 
io _ Yi+2— AYi41 T— 6y; — AY;-1 + Yi-2 
Yi = 74 


cuyos errores de discretización individuales han sido ya expresados. 

Se forma así un esquema en diferencias que relaciona los distintos puntos de [0, [] que serán 
los puntos interiores que deben verificar la ecuación de campo. La ecuación diferencial queda 
aproximada en cada uno de estos nodos de una forma genérica por: 


0% (Rh) Yi-2 + 01 (Rh) Yi-1 +00 (h) y + P1 (A) Yisr + B2[(R) Yiza = fi = $ (x;) 


obteniendo como resultado un conjunto de ecuaciones algebraicas lineales cuyos coeficientes 
dependerán del valor de h y por tanto de la elección de N. A este conjunto de ecuaciones 
hay que añadir las condiciones en los extremos en las que, en caso de aparecer derivadas, las 
aproximaremos utilizando las fórmulas (1.14). 

El conjunto total de estas ecuaciones más las condiciones en los extremos debe permitir la 
obtención de los valores aproximados de y (x,) de la solución buscada en los puntos de |[0, /] 
incluidos, en su caso, los extremos. El sistema de ecuaciones deberá ser, por lo tanto, un sis- 
tema compatible determinado que nos permita la obtención de la solución proyectada sobre 
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la red unidimensional uniforme definida para h. La utilización de esquemas de aproximación 
consistentes y estables nos permitirá que para valores de h progresivamente menores obten- 
gamos sistemas con mayor número de ecuaciones que proporcionan una mejor aproximación 
discreta a la solución deseada. 

El procedimiento se utiliza igualmente en el caso de ecuaciones lineales de coeficientes 
variables, de la forma 


as (1) y" (2) + 03 (2) y” (2) + a7 (2) y" (2) + as (2) y (2) + 00 (2) y (2) =$ (2) 
con aj, f € C([0, 1]). En este caso los coeficientes a (h) y B(h) de cada ecuación algebraica 
no permanecen invariables. 

Por otra parte, hay diferentes formas de definir la proyección sobre la red discreta de la 
función f (x) que aparece en el segundo miembro de la ecuación de campo de manera que no se 
limite a las funciones continuas y que abarque por ejemplo a las funciones integrables. Pueden 
ser consideradas en la ecuación diferencial funciones f (1) segmentariamente continuas en 
[0, /] con un número finito de discontinuidades de salto finito en dicho intervalo. Como ya se 
ha anotado en $ 1.1, para obtener la proyección de la función sobre la red bastará considerar 
los valores promedios de f (1): 


f:=f (u;) = : (F(27) +1 (27) 


Finalmente puede establecerse, con mínimas restricciones, una serie de puntos interiores en 
el intervalo y trabajar de la misma forma obteniendo un sistema de ecuaciones algebraicas con 
un paso de red variable h en cada ecuación; es decir: trabajar sobre redes no uniformes. Como 
se ve, independientemente del grado de aproximación obtenido, el método de aproximación 
por diferencias finitas puede ser aplicado con una gran generalidad. 

Obtenido el esquema en diferencias que aproxima la ecuación diferencial, éste debe ser 
consistente en el sentido de que el error sea un infinitésimo con h. Para nosotros será suficiente 
precisar el orden de aproximación general que será O (h?). No obstante, utilizaremos también 
el término complementario de Lagrange, para realizar algunos ajustes sencillos. Por otro lado, 
la diferencia entre los valores de la solución y la aproximada obtenida por diferencias será 
generalmente distinta en cada uno de los puntos por lo que hay que introducir una medida 
del orden de aproximación logrado. Si esta medida tiende a cero con h la aproximación por 
diferencias convergerá a la solución y el esquema será estable. 

Nuestro interés es el de proporcionar una descripción puramente informativa de los méto- 
dos en diferencias con un carácter esencialmente práctico. El problema que nos va a ocupar es 
del análisis de la flexión simple en barras rectas de inercia constante que tiene peculiaridades 
propias. Con las hipótesis que se describirán, estos problemas tienen garantizada la existencia 
y unicidad de la solución de los problemas de contorno que de ellas se deducen. Asimismo, los 
esquemas que formaremos son consistentes y estables con la convergencia asegurada cuando 
h => 0. 

Obtenidos los valores discretos en los puntos de la red establecida se puede tratar de 
reconstruir de forma aproximada la solución continua lo que pone en relación estos problemas 
con los de interpolación. En los problemas aquí presentados cuando se pretenda presentar 
una gráfica de la solución aproximada nos limitaremos a trazar una poligonal asociada al 
espectro de las soluciones discretas. 


1.4 Ecuaciones de 2% orden con valores de la solución en los extremos 


1.4. Ecuaciones de 22 orden. Caso en el que se 
conocen los valores de la solución en los extremos 
del intervalo 


Se trata de obtener una aproximación de la solución particular en [0,/] de la ecuación di- 
ferencial ordinaria lineal de segundo orden cuyos coeficientes son funciones reales de variable 
real, en general continuas 


a2 (1) y" (1) + ar (2) y (2) + a0 (2) y (2) = f(z), 


con aj, f € C([0,1]) y az(1) 4 0 para todo zx € [0,/], que satisfacen las condiciones de 
extremo y (0) = A, y (1) = B, teniendo garantizada solución única. 

Elegido N obtendremos el paso h = E y la red de puntos R] [0,1]: es decir, 1, = 0 +¿h, 
¿€ ZN[O, N]. Son conocidos los valores de la solución para los índices ¿=0e ¿= ÑN, es decir, 
Yo = Ae yy = B. Teniendo en cuenta las fórmulas (1.14), la ecuación de campo quedará 
aproximada para cada 1; por: 


Ur = 2 FU Yi+1 — Yi-1 
ay(1) , +a(a) 2 + ao lady = f, =$ (x;) - 
h 2h 
Con la notación az (1) = a3;, a1 (1;) = 41;, 0 (1;) = Ao, puede escribirse para cada punto 
de la red interior al intervalo la ecuación, 0<1<N, 


(249; — har)Yi-1 + 2(h?ap; — 2095)Y; + (209; + har)Yir = 2 ds 


y las condiciones y (0) = A, y (1) = B (yo = A, yw = B) deberán determinar un sistema 
compatible determinado de N — 1 ecuaciones. Las matrices de coeficientes de los sistemas 
algebraicos resultantes son matrices banda características de estos sistemas más significativas 
y notorias cuando la red es uniforme. 

El error debido a la discretización puede venir expresado a través de los términos comple- 
mentarios de Lagrange en función de ciertos valores de las derivadas tercera y cuarta en el 
entorno de cada punto de la red: 


1 1 : 
el (guay (2; Sp ah) =P ay” (Ly + 6) : 


con a, 6 € (0, 1). 

Si la ecuación es de coeficientes constantes, los coeficientes de la ecuación algebraica aproxi- 
mada se conservan lo que simplifica la construcción del sistema algebraico. Si existen puntos 
interiores de la red en la que f(x) tenga una discontinuidad de primera especie bastará, 
como se ha dicho, con tomar el valor promedio para evaluar f;. 


Ejemplo 1. Consideremos la solución y (1) en [0,3] con y(0) =0, y (3) =0 de la ecuación 
y" (1) — y (2) = f(x) 
donde 
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Figura 1.2: Función f (x) y discretización del intervalo (0,3) para (a) h= 1 y (b) h =0.5. 


Procedemos a dividir este intervalo en NV partes iguales. Aquí se ha comenzado por es- 
tablecer N = 6 definiendo una red de paso h = 0.5, con cinco puntos interiores (figura 
1.2(a)). 

Para cada uno de los puntos x; € [0,3], la ecuación diferencial de campo quedará aproxi- 
mada por 

Yi+1 — 2Yi + Yir  Yi41 — Yi 
he 2h 
donde f; = 5 ii (27) +f (27). Es decir, 


=f (1:15) 


5Yi-1 — Yi + Mira = fi, 


ecuación de aplicación sistemática para los cinco puntos interiores del intervalo. Teniendo en 
cuenta las condiciones de extremo, nos permiten escribir el sistema de ecuaciones algebraicas: 


Yo —8Yy1+3Y2 = 1, con Yo = 0, 
y -8y+3Y3 = 1, 

5Y2 —8Y3+3Y4 = 0, 

5Y3 —8Ya+3Ys = —l, 

Ya — 8Y5 + 3Y6 = —l, con Yyó = 0, 


que de forma matricial se puede escribir como: 


-8 3.0 0 0 Yi 1 
5 8.30 0 Ya 1 
05 -—88 3=0 uv l.= 0 
0 0=5 -—8 3 Ya —1 
0. 0 0>=5 -—8 Y5 —1 
y cuya solución es: 

Yo 0.0000 

Yi —0.2631 

Ya —0.3684 

ya | = | —0.2105 

Ya 0.0526 

Y5 0.1579 

Y6 0.0000 


Una mejor aproximación la obtendremos disminuyendo el paso h de la red. Así, si estable- 
cemos una red con h = 0.25 y 11 puntos interiores numerados tal y como muestra la figura 
1.2(b), la ecuación (1.15) quedará ahora: 
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18Y-1 — 32Y4; + 14Yi41 = £: 


y el sistema de ecuaciones resultantes, teniendo en cuenta las condiciones de extremos, será: 


18y, — 32y1 + 14y =  1l, con Yo = 0, 
18y, — 32y2 + 14y3 = 1, 
18ya = 3243 = 14ya = 1, 
18y3 = 32Y4 TF 14ys = 1, 
18y4 — 3245 + 14y8 = 1, 
18y5 — 3246 +14y7 = 0, 
18y4 — 32y7 + 14yg = —l, 
18y7 = 3248 =F 14yo = —1, 
1843 — 32Y9 + 14Yy10 = —1l, 
18Y9 — 32y10 + 14yi1 = —l, 
18410 — 32411 + 14y2 = —1, con y12 = 0, 


De forma matricial el sistema se puede escribir ahora como: 


32 14.0 0 0 0 0 0 0 0.0 Y 1 
18 0. 0. 0. 0. 0. 0.0 Ya 1 
0 0. 0. 0. 0. 0. 0.0 Ya 1 
0 0 1 0. 0. 0. 0. 0.0 Ya 1 
0 0 M0 0. 0 0.0 Ys 1 
0. 0.0 0. 0. 0.0 y |=| 0 
F y 6 0 0. 0.0 Ya =i 
0. 0. 0. 0.0 0. 0 Ys =i 
0. 0. 0. 0. 0.0 0 Yo | 
0. 0. 0. 0. 0.0 Y1o =i 
FC DU Ó 0 0 0 0 0 0518 <-32 Yi il 


y la solución del sistema es: 


Yo 0.0000 
Y (0.1464 
Yo 0.2632 
Ya 0.3420 
Ya 0.3719 
ys —0.3389 
ys |= | —0.2250 
yr —0.0786 
Ys 0.0382 
Y10 0.1469 
Yi 0.1139 
Y1 0.0000 


Las poligonales que aproximan la solución buscada en ambas casos pueden representarse 
superpuestas como se muestra en la figura 1.3(a) junto con la solución exacta que en este 
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0.3 7 7 7 ! r 10 
exacta 


021 —3—h=0.5 
—2— h=0.25 


error 
E 
o 


Figura 1.3: (a) Solución analítica junto con la solución numérica para h = 0.5 y h = 0.25; (b) gráfica de error 
junto con gráfica de la función C'h?, con C'= 1. 


caso particular puede calcularse mediante los métodos clásicos de resolución de ecuaciones 
diferenciales. Se aprecia de forma clara como el error de la solución aproximada es menor 
cuanto menor es el valor de h. 

En general podemos escribir el sistema de forma matricial en función del parámetro h 
como: 


2 l=5 Y h 
14h 2 104 0 , 
1 14+5 e L=a Ñ Yy3 f3 
En Ls ma. LE ya = da 
14% -2 1-£ YN-2 Ín-2 
14%  -2 YN-1 Ín-1 


siendo nulos los coeficientes de la matriz que no aparecen representados. 

Si resolvemos el problema para distintos valores de h, y calculamos el error máximo que 
se comete en cada punto de R) [0, l], podemos representar una gráfica de error, figura 1.3(b), 
en la que cada círculo representa el error cometido para cada valor de h. En la figura se 
muestra además la función Ch?, con C = 1 en este caso, como referencia, que en escala 
logarítmica viene representada como una recta paralela a la de la función de error, pues tal 
y como habíamos adelantado el orden de error del método es O (h2). 


1.5. Ecuaciones de 22 orden. Caso de la presencia de 
derivadas en las condiciones de extremo. Puntos 
“ficticios” 


Si en el intervalo [0,/] las condiciones frontera involucran el valor de la derivada y no el de 
la función, este último pasa a ser una incógnita más en el sistema de ecuaciones algebraicas 
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resultante pero en este caso dispondremos de una ecuación adicional en el extremo, lo que 
nos obligará a considerar formalmente valores de la solución fuera del intervalo que en la 
práctica supone la extensión de la red unidimensional R; [0,1] a uno o ambos lados del 
intervalo. Á estos puntos exteriores alcanzados por el esquema en diferencias que aproxima 
la ecuación diferencial se les puede denominar “puntos ficticios”. La ecuación adicional que 
liga los puntos ficticios de la red con los interiores es consecuencia de utilizar en los puntos 
extremos las correspondientes aproximaciones teniendo en cuenta las fórmulas (1.14). 

En el caso de la presencia de derivadas en alguno de los extremos y su substitución por dife- 
rencias, además de los errores de discretización de la aproximación deben de considerarse los 
de aproximación de estas condiciones adicionales. Estos últimos se “propagan” normalmente 
desde los extremos afectados hacia el interior. 

Tomaremos como ejemplo la misma ecuación diferencial del epígrafe anterior con una 
condición en el extremo derecho que involucre la derivada y/ (1). 


Ejemplo 2. Se trata de obtener ahora una solución aproximada y (x) en [0,3] con y (0) = 1, 
y (3) =0 de la ecuación 


donde 
$ 1 sixe[0,15), 
His) = Í =1,. sire(L5,3l, 


Elegido el paso de red h = 0.5, el sistema de ecuaciones algebraicas será ahora: 


Yo —8y1+3Y2 = 1, con Yo =1, 
5y -8y9+3Y3 = 1, 
5Y2 —8Y3+3Y4 = 0, 
5Y3 — 8Y4 +3Y5 = —1, 
Ya — 8Ys +3Ys = —l, 
5Y5 — 8Y6 + 3Y7 = —l, 
con la condición adicional 
Y7 Ys _ 0 
2h ' 


es decir, y7 — Ys = 0, que involucra un punto ficticio 17 fuera del intervalo a su derecha 
(figura 1.4(a)) que puede ser evaluado en función del punto interior más próximo al extremo, 
es decir, en este caso y7 = Ys. El sistema de ecuaciones resultante quedará por lo tanto: 


-8y1 +3Y2 = -—4 
5y1 —8y9+3y3 = 1 
ya —8Ya+3Yy4 = 0 
5Y3 — 8ya+3Ys = —1 
Ya — 8Y5s + 3Y8 = —1 
8Yys — 8Y8 = —1 
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Figura 1.4: Función f (x) y discretización del intervalo (0,3) para (a) h = 1 y (b) h = 0.5, incluyendo un nodo 
ficticio a la derecha del intervalo. 


que puede expresarse de forma matricial como: 


83.0. 00.00 Yi —d4 
5 8 3 0 0-0 Ya 1 
0. 5 8.3 0-0 Ya | 0 
0 0 5 =8. 3 0 ya | | -—1 
0.0. 0. 5 3 3 Y5 =1 
0.0. 0 0 8 -—8 Y6 -1 
y cuya solución es: 

Yo 1.0000 

Yi 0.7588 

Ya 0.6902 

ya | = | 0.9092 

Ya 1.2742 

Y5 1.5492 

Y6 1.6742 


Para una retícula de paso h = 0.25, el sistema de ecuaciones algebraicas será ahora: 


18y, — 32y1 +14y = 1, con Yo = 1 
18yí — 32y2+14y3 = 1, 
18y2 — 32y3 + 14y4 = 1, 
18y3 — 32y4+14y5 = 1, 
18y4 — 32y5 + 14y = 1, 
18y5 — 32y4 + 14y7 = 0, 
18Y6 — 32y7 + 14yg = —1, 
18y7 —32ya + ly = =l, 
18Yy3 — 32Y9 + 14Y1o = —l, 
18Y9 — 32Y10 + 14Y1 = —1, 
18Y10 — 32Y11 + 14y12 = —1, 
18y11 — 32y12 + 14y13 = —1, con Y13 = Y11 


donde ahora el punto ficticio que queda fuera del intervalo será 113. Este sistema de 12 
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ecuaciones expresado de forma matricial será: 


-32 14 0 0 0 0 0 0 0 0 0 0 Yi 17 
18 -—32 14 0 0 0 0 0 0 0 0 0 Ya 1 
0 18 -—32 14 0 0 0 0 0 0 0 0 Y3 1 
0 0 18 -—32 14 0 0 0 0 0 0 0 Ya 1 
0 0 0 18 -—32 14 0 0 0 0 0 0 Y5 1 
0 0 0 0 18 -—32 14 0 0 0 0 0 Y6 0 
0 0 0 0 0 18 -—32 14 0 0 0 0 Y7 o -1 
0 0 0 0 0 0 18 -—32 14 0 0 0 Y8 -1 
0 0 0 0 0 0 0 18 -—32 14 0 0 Yo -1 
0 0 0 0 0 0 0 0 18 -—32 14 0 Y10 -1 
0 0 0 0 0 0 0 0 0 18 -—32 14 Y11 -1 
0 0 0 0 0 0 0 0 0 0 32 -—32 Y12 —1 
y cuya solución es: 

YO 1.0000 

Yi 0.8627 

Yo 0.7576 

Y3 0.6930 

Ya 0.6835 

Y6 = | 0.8876 

Y7 1.0753 

Ys 1.2452 

Yo 1.3923 

Y10 1.5100 

Y11 1.5897 

Y12 1.6211 


Las gráficas de ambas soluciones pueden compararse en la figura 1.5(a). 
Matricialmente podemos escribir el sistema resultante para un valor de » arbitrario como 


AY =F, 
donde 
2 1-£% 
h h 
1+4 2 1-% , 
14% -2 1-4 
1 14% -2 1-2 
A ho _h ; 

h 1+* 2 1-* 
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(a) (b) 


2 , , 7 - : 10 
exacta 


Le l| 9 h=0.5 
—2— h=0.25 


error 
a 


Figura 1.5: (a) Solución analítica junto con la solución numérica para h = 0.5 y h = 0.25; (b) gráfica de error 
junto con gráfica de la función C'h?, con C'= 1. 


Yi l— > (1 + 2) 
Ya 
Ya 1 
Y YN/2-1 ! P- pl 
UN/2 0 
YN/2+1 =i 
YN-1 sal 
YN —1 


La gráfica de errores para distintos valores de h (calculados como en el ejemplo anterior) 
se muestra en la figura 1.5(b). 


1.6. La ecuación y” (1) = f (1) 


Particularmente nos interesa abordar mediante diferencias la integración de la ecuación 
y"(x) = f (xv) en el caso de determinadas funciones f (1) segmentariamente de clase cero 
en J = [0,/], con un número finito de discontinuidades de primera especie que, además de 
los saltos finitos, presentan un número finito de quebraduras. Además, las condiciones adi- 
cionales pueden venir definidas como condiciones de extremo o bien incluir alguna condición 
en algún punto del interior del intervalo. Aparte de las discontinuidades de la propia función 
f(x), los puntos en los que esté definido el ángulo de una línea quebrada, por ejemplo, pro- 
porcionarán también una cierta contribución al error de discretización local y a estos errores 
habrá que añadir los de las condiciones adicionales si estas deben ser asimismo aproximadas. 
En todo caso consideraremos que las anomalías de la función f (1), ya sean discontinuidades 
o quebraduras coinciden con los nodos de la red discreta sobre la que se proyecta. 

En estas condiciones, la aproximación de la ecuación diferencial resulta sencilla. Para cada 
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(b) 


x=3 
(12345678€9101 


M2 13 14 15 16 17 12 19 
x=9 


Figura 1.6: Función f (x) y discretización del intervalo (0,9) para (a) h = 1 y (b) h = 0.5, incluyendo un nodo 
ficticio a la derecha del intervalo. 


punto x; de le red unidimensional que cubre el intervalo se deberá verificar simplemente: 


Visa 244 Yi 
J3 


= fi 


y el error de discretización local puede evaluarse por el comportamiento del factor 
5 1y" (a; +0h) con 0 < 0 < 1, factor que depende del comportamiento de y” (%) = f” (x) 
en el entorno de cada nodo zx; por lo que habrá que abordar la función f (x) a la derecha o 
la izquierda para evitar o bien sus discontinuidades o bien las de su derivada. 


Ejemplo 3. Se trata de encontrar la solución particular y (1) de clase uno, que en el intervalo 
[0, 9] satisface la ecuación diferencial y” (1) = f (1), donde 


SS —z, si x E [0, 3) 
ol 2 si a € [3,9] 


con las condiciones adicionales y (3) =0 e y (9) =0. 

La función f (1) es una función continua con una quebradura en el punto x = 3 en el 
que su primera derivada es discontinua. La solución buscada es desconocida en el extremo 
izquierdo del intervalo pero es conocida justamente en el punto x = 3. 


Si consideramos, por ejemplo, una red de paso h = 1, figura 1.6(a), aplicaremos en los 
puntos interiores la ecuación 


Yi-1 + 2Yi + Yira = fi, 


pro ye de la que se deduce que Y¡o = Ys- 


2h 


con la condición en el extremo derecho 0 = yy = 
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Se obtendrá, por lo tanto, el sistema: 


Yo =2Y +ya = —L, 

Y =2Y+Y = =2, conyg=0 — Y 2y2= -2, 
y —2Y3+Y. = -—3, conyg=0 —> YaT+Y=-3, 
Ya —2Ya+Ys = -—2, conyg=0 —> —2Ya4 + Ys = -—2, 
Y —2Ys+Y = —L, 

Ys —2Y6+Yy7 = 0, 

Ye —2y7+Y8 = 1, 

yr —2Yg+Y = 2, 

Ys —2Y9+Yo = 3, con Yo =Yg8 =>  2Yg — 2Y9 =3, 


sistema compatible determinado que puede expresarse matricialmente como: 


1-2 1 0 0 0 0 0-0 YO —1 
0 1 20. 0 0 0 0.0 Yi —2 
0-0 1 1 0 0 0 0 0 Ya —3 
0.0 0 -—2 1 0. 0 0-0 Ya —2 
0.0 0 1 -2 1 0 0 0 y l|=| -1 
0.00 0.0 1-2 1 0 0 Y6 0 
0.00 0 0-0 1-2 1 0 Y7 1 
0.00. 0 0 0.0 1-2 1 Y8 2 
0.0. 0 0 0 0 0 2-2 Yo 3 
con solución: 

YO —9.5000 

Yi —5.0000 

Ya —1.5000 

Y3 0.0000 

Ya | —1.5000 

ys | —5.0000 

Y6 —9.5000 

Y7 —14.0000 

Y8 —17.5000 

Yo —19.0000 


Una red más estrecha sobre el intervalo, figura 1.6(b), nos permite ajustar un poco más 
el conjunto de soluciones. Utilizando una discretización de paso h = 0.5, por ejemplo, el 
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) 
) 
) 
) 
) 


Yo — 2Y1 + Ya 
Y1 — 2Y2 + Y3 
Ya — 2Y3 + Ya 
Y3 — 2Ya + Ys 
Ya — 245 + Yo 


sl 
zl 
zl 
zl 
sa 
e 
pl 


E > E Es 3 
a la lo lo lo 


sistema de 18 ecuaciones es: 


ja + 
mia | 
! 
2 SS 
D» 
F+ | 
iS E 
DMA D 
AN 


= 
= 
= 


=.0 
=.0 
Yi7 


con Y19 


) 
) 
) 
) 
) 
) 
) 
) 


SA A A AA 


Laa... 
DADAS 
ÓN 
o . 0 o TN mios ano Nros 
DZDAIDAID ADD 
ANA A AAA NANO 
A O A O A 
le 
Ss 3355558585 
Ep DDD DXDDXOD DDD 


| 5] ¡SN a] ¡SN 5) 5) 5] 
2 Epi Sa ll pl ol pl e El 
o 3 la la la la la la 


que matricialmente viene expresado mediante: 
AY =F 


donde 


oc0o00o00o00o0000000000000 rr 
sosososnosoosasoaa: 


oco0o0o00o00o00o00o000000000 py A 


NnO00O00O000O00000o000o0o00o00nsN> 
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Yo —0.5 
Y1 =1 
Ya 1.5 
Y3 =2 
Ya —2.5 
Y5 =3 
Y7 —2.5 
Ys =2 
Yo 1.5 
ds Y10 il ás ==] 
Yi1 0.5 
Y12 0 
Y13 0.5 
Yia 1 
Y15 1.5 
Y16 2 
Yi7 2.5 
Y1i8 3 
y que tiene como solución: 

YO 9.1250 

Yi —6.8750 

Ya —4.7500 

Y3 2.8750 

Ya —1.3750 

Y5 —0.3750 

Y6 0.0000 

Y7 —0.3750 

Y8 —1.3750 

Yo ES —2.8750 

Y10 —4.7500 

Y11 —6.8750 

Y12 9.1250 

Y13 —11.3750 

Ya -13.5000 

Yis —-18.2500 


La solución exacta de esta ecuación puede ser obtenida sin ninguna dificultad por métodos 
continuos viniendo definida por segmentos. En la tabla 1.1 mostramos una comparación entre 
la solución exacta y la aproximada obtenida para h=1yh=0.5: 


Igual que en ejemplos anteriores, podemos escribir el sistema a resolver para un valor de 
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1.6 La ecuación y” (1) = f (x) 


tel R=L | R=053 ute 

0.0 | -9.5000 | -9.1250 9.0000 
0.5 -6.8750 6.7708 
1.0 -5.0000 | -4.7500 | 4.6667 
1.5 -2.8750 | 2.8125 
2.0 | -1.5000 | -1.3750 1.3333 
2.5 -0.3750 |. 0.3941 
3.0 0.0000 0.0000 | 0.0000 
3.0 -0.3750 | 0.3941 
4.0 | -1.5000 | -1.3750 | 1.3333 
4.5 -2.8750 | 2.8125 
9.0 | -5.0000 | -4.7500 4.6667 
9.0 -6.8750 6.7708 
6.0 | -9.5000 | -9.1250 9.0000 
6.5 -11.3750 | 11.2291 
7.0 | -14.0000 | -13.5000  13.3333 
7.5 -15.3750 | 15.1875 
8.0 | -17.5000 | -16.8750 | 16.6667 
8.5 -17.8750 | 17.6458 
0.0 | -19.0000 | -18.2500 | 18.0000 


Tabla 1.1: Comparación entre la solución exacta y la aproximada obtenida para h= 1 y h= 0.5. 


h genérico, en este caso, el sistema matricial a resolver será AY = F', con: 


1-2 1 
1 -—2 1 
1 -—2 1 
1 -—2 
1 
7 1 1 
2 ) 
h =2 1 
1 2 1 
1-2 1 
2 -—2 
Yo f (20) 
T 
Y — a Si He Fes) 
UN f (xn) 


En este caso y puesto que la función f(x) está formada por rectas cuya pendiente es 
constante, los errores se introducen localmente en la quebradura sobre el punto x= 3 y en 
el extremo derecho por la aproximación de la condición frontera (es decir, en las ecuaciones 
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3 x=3 
M N 0607 08 09 


15 16 17 18 19 
x=9 


Figura 1.7: Puntos ficticios añadidos para mejorar el error de discretización: (a) a la derecha del intervalo y (b) 
en 2=3. 


6% y 18%, de extremo, para h = 0.5), toda vez que y*” (x) en el resto de los puntos es 
nula a su derecha e izquierda. Con el fin de ilustrar el comportamiento de los errores y los 
posibles ajustes y aproximaciones, podemos corregir en este caso los errores cambiando las 
dos ecuaciones que introducen los errores locales. 

Con respecto a la condición del extremo derecho es a la izquierda de este punto yg (1) 
que tiene sentido físico la ecuación general y la condición adicional de extremo, esto es: la 
ecuación a aplicar es correctamente (utilizando desarrollos de Taylor en 117 y 119): 


O mE 0) a (1.16) 


en tanto la condición de extremo se escribirá: 


epa Ys—=Yr 11, 
Y1g (1) = 2-05 6/0 


y puesto que y¡g (¿) es constante e igual a 1 determinará y19 = Y17 + 3, (figura 1.7(a)). Este 
valor llevado a (1.16) nos permite escribir la ecuación: 


17 
2Y17 — 2Y1i8 = 55» 


24 
que es la ecuación que debe de sustituir a la 18* en el sistema. 

Si la solución es de clase uno, tenemos garantizada la igualdad de la derivada primera a 
derecha e izquierda de los puntos interiores. El procedimiento es aplicable para el ajuste de 
la ecuación 6%. Abordaremos el punto x = 3 por su derecha e izquierda con puntos (ficticios) 
M y N respectivamente (figura 1.7(b)). 

Por la izquierda tendremos: 


AS 04 
y además 
YM + Y5 
A 
0,5% 
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(b) 


—S— error 
=—— Ch? 


4 3 2 A 0 


10 10 10 10 10 


Figura 1.8: (a) Solución analítica junto con la solución numérica para h = 1 y h = 0.5; (b) gráfica de error junto 


con gráfica de la función C'h?, con C = 10 


de donde 


Por la derecha: 


0D 64d. 
y además 
YI FYN __3 
Mas á 
de donde 
¿(d) = 2y7 + sl 
Y6 = 447 94 
y la igualdad de ambas derivadas: 
17 17 
== 2 pe! 
Y5—= 27 Y + 24 
nos proporciona la ecuación 
17 
2 2Y7 == 
Ys + 2Y7 12 


que es la ecuación que debe de sustituir a la 6* del sistema. 

Las soluciones del nuevo sistema de ecuaciones corregidas la 6% y 18% coinciden ahora con 
las exactas. 

En general este tipo de ajustes no se realizan o bien directamente no pueden realizarse y 
se prefiere obtener una mayor aproximación reduciendo el paso h de la red. 

En la figura 1.8(a) se representan tres gráficas: la poligonal obtenida para la aproximación 
sobre red h = 1, la obtenida para la aproximación h = 0.5 y la solución continua exacta, 
mientras que en la figura 1.8(b) se muestra la gráfica de errores obtenida al resolver el 
problema para distintos valores de h. 
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1.7. Ecuaciones de 4% orden. La ecuación y” (1) = f (x) 


Establecida una discretización de paso h que cubre el intervalo J = [0,[], la substitución 
de las derivadas por sus aproximaciones en diferencias centrales O (h?) en la ecuación: 


a4 (2) y” (2) + as (2) y” (2) + ar (2) y” (2) +01 (2) y (2) +00 (2) y (2) =f (2), 


con as € C%(J), j = 0,1,2,3,4, dará como resultado una ecuación algebraica para cada 
punto interior 2;: 


be Yi+2 — AYi41 + 6Yi — ÁYi-1 + Yi-2 
4i 


Po — Zar + 2Yia1 — Yi-2 
ht y 


2/8 


VYi41— 25 + Yi E VYi41 — Yi-1 
21 p2 I 1: 2h, 


Se deberá disponer, además, de cuatro condiciones adicionales que garanticen la existencia 
de una solución única del problema. En general, aparecen un par de condiciones en cada 
extremo, expresadas por ecuaciones que involucran los valores de la solución y de algunas 
de sus derivadas. Nuestro interés se centra en encontrar soluciones en [0, /] de las ecuaciones 
más sencillas del tipo y” (1) = f (x) donde en general la función f (1) tiene, en el intervalo, 
un número finito de quebraduras o de discontinuidades de primera especie. 

En el caso particular que tratamos aquí, discretizada la variable, cada punto del interior 
deberá satisfacer la ecuación de campo aproximada: 


Eduti= $e 


Yi+2 — AYiz1 + 6Yí — AYi-1 + Yi-2 
pa 


AS HER 
con 


f=>3 (4 (07) +1 (07) 


y donde las condiciones en los extremos involucran las derivadas de la función buscada por lo 
que, independientemente de ser conocidos en la frontera los valores de la solución, existirán 
puntos exteriores alcanzados por el esquema en diferencias o “puntos ficticios”. 


Ejemplo 4. Se trata de obtener la solución que satisface en el intervalo [0,4] la ecuación 


q (e) =F (0), com: a 
d, six € |0,2), 
fla, sia € (2,4), 


y las condiciones en los extremos: y” (0) =1, y” (0) =1, y (4) =0, y' (4) =0. 
Comenzaremos por definir una discretización de paso h. Elegimos, por ejemplo, la de paso 
h = 0.5 y numeramos los puntos de la red en el interior tal y como se representa en la figura 


1.9(a). 
Las condiciones en el extremo izquierdo podrán aproximarse por: 


mo Ya — 2Y1 =P 2Y-1 —= Y-2 e: 
Yo — 2h3 — 
Yi — 2Yo + Y-1 
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(b) 


al iii) PORN 161 
+ 


x=2 x=2 


Figura 1.9: Función f (x) y discretización del intervalo (0, 4) para (a) h = 0.5 y (b) h = 0.25, incluyendo un nodo 
ficticio a la derecha del intervalo. 


que nos obligan a extender por la izquierda los puntos ficticios 1_1 y 1_2 fuera del intervalo. 
Los valores de la solución en estos puntos exteriores quedan relacionados con los interiores 
resolviendo las dos ecuaciones y obteniendo: 


1 
Ya = Yo Y += Yo Y + y 


il 
Ya = 2y-1— 21 + Ya — 2h? = Ayo — dy + Ya +2h% — 2h* = dyo — dy + Yo + y 


y en el extremo derecho, además de yg = 0, la condición adicional 


/ Yo — Y7 
+ E O, 
Yg 2h 
permite asignar al punto 19 “ficticio” el valor yg = yz. 


El sistema de ecuaciones para h = 0.5 quedará por lo tanto: 


1 y-2 = Ayo — 4Yy1 + Ya +3, 
a (Y-2 — 4Yy-1 + 6Yo — 4y1 + = 0, con Ñ 
057 (Y-2 — 4Yy-1 + 6yo — 4Yy1 + Ya) l ya=2Y0-4+1, 

03 (Y-1 — Ayo + 6Y—4ya+ Ys) = 0.5, con Yy1=?2Y0-Y+; 

a (Yo — M1 + 6y2 — Aya + ya) = 1, 


== (Yy1 — 4y2 + ya — 4Ya + Ys 


( ) 

05 (Ya — 4y3 + 6ya — 4ys + Yo) 

03 (Ys — Aa + 6y5 —4Ye+Yy1r) = —1.5, 
( ) 
( ) 


4 


E (Ya — As + 6Y6 — Ayr + Ys) = —1, con y¿=0 


05 
= (Ys — 4yó + 6y7 — ys + Yo) = —0.5, con ys=0, yo =yr 


0.5% 
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que de forma matricial se puede escribir: 


2 -42 0 0 0 0 0 Yo 0 1-4 
2.5 4 10 0 0.0 Yi 0.5 5 
1-46 4 1 0 0 0 Ya 1 0 
1 0 1 -4 6 -4 1 0 0 Ya | 1.5 il 0 
0o5t| 0 0 1 -4 6 -4 1 0 yal | 0 | 05t| 0 
0. 0 0 1 -4 6 -4 1 Y5 —1.5 0 
0.0 0 0 1 -4 6 -4 Y6 —1 0 
0.0 0 0 0 1-4 7 Y7 —0.5 0 
y cuya solución es: 
Yo 44.1250 
Yi 35.7500 
Yo 27.7500 
Y3 20.2812 
ya | =| 13.5625 
Y5 7.9062 
Y6 3.6250 
Y7 0.9375 
Ys 0.0000 


Estableciendo una red mas fina, con h = 0.25, como se representa en la figura 1.9(b), el 
sistema será ahora: 


3 
sa [Y-2 — 4y-1 + 6Yo — 4Y1 + Ya) = 0, con Í $ _ e Ñ Ao + 32 
027 3 (Yi — Yo +6y—4y+ ys) = 0.25, con Yy1=2Y0-Y+>5% 
025 (Yo — 4y1 + 6ya — 4yg+ ya) = 0.5, 
E 7 (y — Aya + 6y3 — 4ya + Ys) = 0.75, 
3% 3 (Yo — 4Y3 + 6Ya — 4Y5 + Y6)  = E 
O. 57 (y3 — 4ya + 6ys — 4y6 + yz) = 1.25, 
0. 25 (Ya — 4ys + 6yé — 4y7 + ys) = 1.5, 
025% (Ys — 4Y6 + 6y7 — 4yg + Yo)  = 1.75, 
025 (Ys — 4yr +6ys — y + Yo) = 0, 
035 (yz — 4Ys + 6Yo — 4Yio + Ya) = —1.75, 
035 (Ys — 4Yo + 6Y1o — AY +Yia) = —1.5, 
05 (Y — Ayio + 6ya1 — Aia + Yis) = —1.25, 
ba (Yo — Ayu + 6y12 — Ayiz + Ya) =  —1l, 
027 (ya — 4y12 + 6yi3 — Ayra + Yas) = —0.75, 
sd (Y12 — AYi3 + 6yra — 4Yi5s + Y16) = 0.5, con Yyió6=0 
03% 2 (Y1z — 4Yia + 6Yis — AYió6 + Yir) = —0.25, con Yi6=0, Yir = Ys 
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con: 


220200000000 O0ri y 


220000000000 r ¡y Y 


E 
=|8 

Y 
¡Booocooocnocoooo 
| 
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y cuya solución es: 


Yo 44.6328 
Yi 40.3516 
Ya 36.1484 
Y3 32.0400 
Ya 28.0449 
Y5 24.1846 
Y6 20.4844 
Y7 16.9746 
Ys =. | :138:6914 
Y9 10.6777 
Y10 7.9766 
Y11 5.6240 
Y12 3.6504 
Y13 2.0810 
Y14 0.9375 
Yi5 0.2383 
Y16 0.0000 


En general, podemos escribir el sistema para un valor de h genérico de la forma: 


AY =F, 
donde en este caso, 
2 -4 2 
2.5 -—4 1 


h4 = Ñ ») 
1-4 6 -4 1 
ll -4 6 -4 
1 -4 7 

Yo f (20) —2h* — 2h* 
Yi fx1) h? 
Y2 f (22) 0 
Y3 f (x3) 1 0 
Y= : PS . = hs ] 
YN-=3 Fltn-3) 0 
YN-2 E (1-2) 0 
YN-1 6d (2-1) 0 


Las gráficas de ambas aproximaciones pueden compararse en la figura 1.10(a) mientras 
que en la figura 1.10(b) se muestra la gráfica de errores obtenidos para distintos valores de 


h. 
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(a) (b) 


50 : : 10' 
exacta —S— error 
: —8— h=0.5 o [| ———cCh? 
ao 0 —2— h=0.25 |. 15 
107 
30 
10% 
20 
107 
10 10 
0 107 3 2 1 0 
0 107 107 107 10 


Figura 1.10: (a) Solución analítica junto con la solución numérica para h = 0.5 y h = 0.25; (b) gráfica de error 
junto con gráfica de la función Ch? con C=1. 
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